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1) Introduction 

Very soon after the formation of a neutron star, thermal energies inside it 
drop below the energy gap associated with neutron pairing. One therefore 
expects all but the outermost parts of the crust to pass rapidly into a su¬ 
perfluid state il, 0- The critical angular velocity for vortex formation is 
negligibly small for bodies of astronomical size and therefore, as hrst noted 
by Ginzburg and Kirzhnits [Q, the star’s interior will be threaded by vortex 
lines parallel to the axis of rotation. 

A stimulating review by Kirzhnits and Yudin discusses some of the in¬ 
teresting and surprising general-relativistic effects associated with superfluid 
neutron stars. Consider, for instance, an imaginary experiment in which the 
star is initially at rest and a spin is imparted to the crust. Tense-Thirring 
(frame-dragging) effects will make it appear that the superfluid interior has 
been set into partial co-rotation, i. e. , the contravariant azimuthal veloc¬ 
ity vf = d(l)/dt with respect to distant stationary observers will become 
non-zero. If the crust is separated from the fluid by a small gap, then this 
induced rotation can easily exceed the critical angular velocity without form¬ 
ing any vortices. This is because the superfluid’s orbital angular momentum, 
proportional to the covariant component Vg^,, remains zero. The induced ro¬ 
tation also makes no contribution to the value of the metric coefficient g^t 
(the “gravimagnetic potential”), whose non-vanishing is due entirely to the 
rotation of the crust. 

Superfluid flow is in all circumstances irrotational. For regular axisym- 
metric flows this implies Vgrp = 0, i. e. , the superfluid can never acquire 
orbital angular momentum. It is nevertheless able to entrap angular momen¬ 
tum within the cores of its vortex hlaments, which consist of normal fluid. 
To simplify the description, one generally assumes that it is permissible to 
replace the vortex array by a hctitious “effective” orbital flow ^ 0) 

when averaging over macroscopic domains, i. e. , domains large compared 
with the vortex spacing (about 10“^cm for the Crab pulsar). The equiva¬ 
lence assumption is that gravitational and other physical effects of the vortex 
array are macroscopically indistinguishable from an effective orbital motion 
having the same distribution of angular momentum. This seemingly inno¬ 
cent assumption is called into question in a recent analysis by Kirzhnits and 
Yudin P . The authors do not succeed in hnding an acceptable solution rep¬ 
resenting the contribution of a vortex line to the external gravitational held 
of a spherical body in the linearized Einstein theory. They conclude: “The 
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general formula relating the angular momentum of a body with the asymp¬ 
totic behavior of its gravimagnetic field is invalid when vortex hlaments are 
present.” 

This conclusion is not easy to understand intuitively, since the formation 
and slow migration of vortices into the superfluid from the interface with the 
crust is macroscopically an axisymmetric process, which conserves angular 
momentum locally and involves no appreciable redistribution of mass. It is 
therefore hard to see how it can affect the external gravitational held at all, 
much less suppress all gravimagnetic effects, as the angular momentum of the 
crust gradually becomes “vorticized”. In short, there is no obvious reason 
why spin and orbital angular momentum should behave in a fundamentally 
different way as sources of gravity. 

The issue is perhaps elementary but clearly important because of the po¬ 
tential implictaions for the gravitational dynamics of the binary pulsar and 
gravitational radiation from neutron star mergers. To remove any possible 
doubt it seems worth devoting some effort to an explicit verihcation that 
bulk rotation and vortex hlaments are indeed equivalent in their gravita¬ 
tional ehects, and this is the object of our paper. In section 2 we review the 
phenomenological form taken by the Einstein held equations for a medium 
which is gravitationally polarized by alignment of internal spins. The ehect 
of this spin polarization on the gravimagnetic potential near inhnity can be 
inferred from an integral identity due to Komar, and this is discussed in sec¬ 
tion 3. These results are exact and valid for strong helds. If the gravitational 
held is weak (i. e. non-relativistic) the linearized Einstein equations provide a 
complete solution of the problem (section 4). Finally section 5 is devoted to 
the simplest realisation of the superhuid neutron star within the framework 
of linearized gravity, namely a spherical star of constant density with purely 
azimuthal, locally irrotational how. This how corresponds to having a single 
vortex hlament along the axis. The origin of the contribution to the far-held 
gravimagnetic held is traced back to the singularity of the how along the axis 
of rotation. 

2) Einstein field equations for a medium with 
internal spin 

Analogously to the Maxwell-Lorentz equations for a dielectric medium, there 
is a statistically averaged form of the Einstein held equations appropriate for 
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the coarse-grained description of a gravitationally polarized granular medium 
1^. In the case of interest to us here the polarization is a spin alignment of 
the granules. 

At the sub-macrolevel, the material inside each granule is described by 
a symmetric conserved tensor representing the flux of 4-momentum, 

which couples to the Einstein tensor in the standard way. On larger scales, 
it becomes meaningful to split the complex motion of the material into a 
translational motion of the granules and a distribution of internal spins. 
Correspondingly the coarse-grained average =< > splits into T“^, 

the flux of translational 4-momentum and a contribution from the flux 
of spin angular momentum (skew in the hrst two indices: 

Detailed analysis shows that the coarse-grained form of the Einstein field 
equations is 


J.ab ^y^jjabc^ jjabc 


1 

2 


^gacb ^ g 


bca 



The translational momentum flux is neither symmetric nor in general 
(in a curved spacetime) conserved. The symmetry and conservation of 
are equivalent to the conditions 

( 2 ) 

Equations (|^ are formally identical with the well-known Belinfante-Rosen- 
feld prescription for symmetrizing what is called (in a field-theoretic context) 
the “canonical” stress-energy tensor [p. 

The simplest illustrative example is a gas of spinning particles. If the skew 
tensor denotes the angular momentum of a particle, its 4-momentum 
and its 4-velocity, then 

T“''(a:) = j N{x,p,s)p‘^v^ujg{p,s), S'“''‘'(a;) = j N{x,p, s)s‘^%''ujg{p, s), 

where N = N{x,p,s) is a distribution function and ujg{p,s) the invariant 
volume form on spin-momentum space. Here the asymmetry of arises 
from the non-collinearity of and for a spinning particle. 
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3) Tolman-Komar integral identities 


To check that the spin angular momentum concealed in (|l|) makes its proper 
contribution to the gravimagnetic potential g^ti h is simplest to appeal to the 
Tolman-Komar integral identities (see e. g. |^. These are rigorous general- 


relativistic analogues of Gauss’s flux theorem, which hold for spacetimes 
admitting Killing-symmetries and relate volume integrals over the matter 
content of an isolated system to surface integrals at inhnity 


V(a^b) =0 ^ 

J = 

as 


= -iaR\ 

-Stt j UT\ - 


( 3 ) 


For stationary helds = df) we have Tolman’s theorem 


where the integration is over at = const 3-space, i runs over the three spatial 
coordinates, and the gravitational mass m is dehned by the asymptotic form 

gtt ^ -(1-) (r ^ cx)). 

r 

For an axisymmetric system, the gravitating angular momentum ma, dehned 
by the asymptotic form 

2ma sin^ 9 . , 

g^pt ~- [r oo) 

r 

in asymptotically spherical coordinates, is similarly given by 


ma 


T/V^d^x, 


( 4 ) 


For our purpose is most conveniently recast in the covariant form 
where is the axial Killing vector. Making use of (|^), the identity 


jjabc _ jjcba _ 


cba 


'^abc 


( 5 ) 


4 





and the Killing property that Va^b is skew, we hnd 

J - rj‘)d^ab = J ^ y 

ESS 

Since the total derivative integrates to zero, (H) becomes 

ma = j j cQV^d^x, ( 6 ) 

in which the separate contribntions of the orbital and spin angular momen¬ 
tum are displayed explicitly. To check that the second term has indeed the 
right form, we go to the non-relativistic weak-field limit, in which we can set 

P 

in cylindrical coordinates is the three-dimensional permutation symbol). 
Then the second term in (H) reduces to / Szd^x, giving the ^-component of 
the total spin angular momentum as it should. 


4) Linearized theory 


In the case of weak helds we can obtain a complete solution to the problem 
by integrating the linearized Einstein equations . 


In approximately rectangular coordinates the metric perturbation hab and 
its trace-reversed form hab are “small” quantities dehned by 

1 , 


9ab hab T h, 


ab 


hnh - h, 


ab 


'-hrjab, 


where rjab is the Minkowski metric. The coordinate freedom allows us to 
impose the four (“harmonic”) gauge conditions dah°’^ = 0, with respect to 
which the Einstein equations and their integrability condition (the linearized 
conservation laws) reduce to the simple form 


Gab - ^d^Kb - STlTab, daT-^ - 0, 


where denotes the Minkowski wave operator. We specialize to stationary 
helds where the solution is 


hab{x) 



(7) 
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and the conservation laws imply that (to linear order) the sonrces mnst obey 
the constraints 


/ 


Jid^x 





( 8 ) 


where J* = T** is the momentnm density and the indices i, k rnn from 1 to 
3. (For instance the last of these follows from the integration of di{x^x^J^) 
over a bonnded sonrce and nsing 9-7 = 0.) 

The external held of a compact sonrce, at distances large compared to its 
dimensions, is thus given by 

K, = 0 , hu = —. 

r 

To obtain the gravimagnetic vector-potential Ai = hu = hu we expand 
\x — x'\~^ in ([^ up to the hrst order in x' and note (|^). This yields Ai in 
terms of the source’s angular momentum Lp 


Ai{x) = 2r ^(x X L)i, U 


J{x X J{x))id^x. 


(9) 


For an axisymmetric source it is convenient to switch to spherical polar co¬ 
ordinates, for which 


Aidx^ = —2r ^(Vx'^eikidx'') = -sin^ 6*70, 

r 

where we introduced the conventional notation = ma. Putting all this 
together, we arrive at the standard form of the linearized stationary exterior 
metric 

ds^ = (1 H - ){dx^ + dy^ + dz^) - sirA 9d(j)dt — (1 - )dt^. 

^ f' 'f 

To exhibit directly the separate contributions of spin and orbital angular 
momentum to the gravimagnetic potential in (^, let us assume the orbital 
4-velocity non-relativistic, so that we can set ^ 1 and neglect To 
this order we may set 


gabc ^ gab^^ ^ 
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which implies S'** = S^^Uk- From (Q) we then hnd to hrst order in n* 

^ g^k ^ Q^^Uk ^ ]_^Q ^ gy_ 

The total momentum current T** in (|l|) therefore decomposes as 

f = pu^+ ]^{d X S)\ 

Substituting in (H) hnally yields the expected result 

L* = J {x X puyd^x + j S^d^x, 

verifying that the angular momentum in vortices contributes in the same way 
as orbital angular momentum to the gravimagnetic potential. 


5) A simple example: gravimagnetic field of 
a single vortex filament on the axis of a ho¬ 
mogenous star 

In order to give a concrete realization of the situation discussed in the pre¬ 
vious sections let us consider the simplest possible model by taking the su¬ 
perfluid velocity to be the irrotational azimuthal flow dip and assuming the 
mass-density to be constant within the star. The equation for the gravimag¬ 
netic potential (within the linearized theory) become 

*d*dA = *d*F = j = 9{a — r)d(p (10) 

The fact that the “velocity” dcp becomes singular at the axis may be conve¬ 
niently described by taking the distributional identity 

dcp= ^dy~ydx ^_^ 2 ^ 

x^ + y^ 

into account, where denotes the tensor-distribution on IR^, corresponding 
to d(j). It is not closed and therefore certainly not exact. This fact will 
actually be important in the description of a vortex hlament along the axis 
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of rotation. Integration of the equation for the gravimagnetic potential eg) 
yields 


d*F = *j = 9{a — r)dr -- 

sin 0 

d9 

*F = ujo + 0(a — r)(r — a) - - dwo = 0 ^ ujq = df 

smO 

F = *df — 6{a — r)(r — a)drd(j) (12) 

where ujq denotes an “integration constant” (solution of the homogeneous 
equation). Since the second term in (|^ naively appears to be closed, one 
might be tempted to discard ujq. However, in the light of (^Tj) we hnd 

dF = 0 = d*df + 9{a — r)(r — a)dr 27i5^'^\x)d‘^x, 

A/ = *d*df = —2'Ji9{a — |;2|)(|;2| — a)-^6^‘^\x). (13) 

This shows that we cannot take ujq to be zero, which would have produced 
a gravimagentic held of compact support, namely within the star. It is pre¬ 
cisely the existence of the hlament, represented by concentrated terms in 
0 ). which gives rise to contributions which are not of compact support and 
therefore maybe detected in the far-held. Since our model is sufficiently ideal¬ 
ized it is possible to solve the Laplace-equation for / exactly. The calculation 
is considerably facilitated by taking a 2 ;-derivative of A/, which simplihes the 
inhomogeneity 


Adzf = 4:Tra5^^\x) — 2Ti9{a — |2;|)(5*'^^(a;), 
and gives rise to the solution 


/=-alog|- + J^ + l 


1 


- V(^ - a) + P -{z + a) log ( —^ + J ——— 
+ \/(z + a)2 -K 


( 14 ) 















Inserting ujq = df into © turns the expression for F into 


F = 




dz - -xdx 

P 


2 y z + a + ^J{z + ay + 
xdx 


+ (^^/{z+~ay~+~^ - \/{z - ay + j - d{a - r)(r - a)drd(j). 


(15) 


Integration of ([l^ gives rise to the potential A 


, z + a + \/(z + ay + z + a n. - -75 -^ 

-4=-!^log- ^ + 

4 2 ; — a + y (^ — 0 )"* + p"* 4 


z — a 


4 - a )2 + — ar — r 


which is not zero in the asymptotic region. 


Conclusion 

In this paper we reviewed the influence of spin and orbital angular momentum 
of a material system on its gravitational held. Within the general framework 
of gravitationally “polarized” matter one can show that both of them con¬ 
tribute to the angular momentum detected by an asymptotic observer. This 
investigation settles the recently discussed issue about the imprint of vortex 
lines found in a superhuid neutron star on the asymptotic gravitational held. 
Specihcally in the linearized approximation of general relativity, one can ob¬ 
tain an explicit solution which clearly exhibits the ehect of the vortex line 
on the gravimagnetic potential. 
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Bn equation (17) of this reference, the factor 1/2 on the right-hand side should be 1/4. 
This has consequences for equations (18) and (19). We thank Yosef Verbin pointing out 
this error. 
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